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ABSTRACT 

Although some of the properties of the caustics in planetary microlensing have been known, our understand- 
ing of them is mostly from scattered information based on numerical approaches. In this paper, we conduct 
a comprehensive and analytic analysis of the properties of the planetary caustics, which are one of the two 
sets of caustics in planetary microlensing, those located away from the central star. Under the perturbative 
approximation, we derive analytic expressions for the location, size, and shape of the planetary caustic as a 
function of the star-planet separation and the planet/star mass ratio. Based on these expressions combined with 
those for the central caustic, which is the other set of caustics located close to the central star, we compare the 
similarities and differences between the planetary and central caustics. We also present the expressions for the 
size ratio between the two types of caustics and for the condition of the merging of the two types of caustics. 
These analytic expressions will be useful in understanding the dependence of the planetary lensing behavior on 
the planet parameters and thus in interpreting the planetary lensing signals. 
Subject headings: planetary systems - planets and satellites: general - gravitatinal lensing 



1. INTRODUCTION 

Microlensing is one of the most powerful methods t hat can 
be used to search for extr asolar planets ( Mao & Paczvnski 
ll991tlGould & Loebl 19921) . Recently, two robust microlens- 
ing de tecti ons of exoplanets w ere reported by iBond et ail 
( 2004 and lUdalski et alJ J2005h . 

The signal of a planetary companion to microlens stars is a 
short-duration perturbation to the smooth standard light curve 
of the primary-induced lensing event occurring on a back- 
ground source star. The planetary perturbation occurs when 
the source star passes close to the caustic. The caustic rep- 
resents the set of source positions at which the magnification 
of a point source becomes infinite. Studies of the properties 
of the caustic are important because the characteristics of the 
planetary perturbations in the microlensing light curve depend 
critically on the properties of the caustic. For example, the lo- 
cation of the perturbation on the lensing light curve depends 
on the location of the caustic. The duration of the perturba- 
tion and the probability of detecting the perturbation are pro- 
portional to the caustic size. In addition, the pattern of the 
perturbation is closely related to the caustic shape. Therefore, 
it is essential to understand the properties of caustics for the 
interpretation of the planetary lensing signals. 

Although some of the properties of the caustics in plane- 
tary microlensing have been known, our knowledges of them 
are mostly from scattered information based on numerical ap- 
proaches. The problem of the numerical approach is that the 
dependence of the planetary lensing behavior on the planet 
parameters of the star-planet separation s (normalized by the 
Einstein ring radius Q%) and the planet/star mass ratio q is 
not clear. There have been several attempts to overcome this 
ambiguity using analytic methods. By treat i ng the pla net- 
induce d deviation as a perturbation, Dominik ( 1999) and I Anl 
( 120051) derived analytic expressions for the locations of the 
central caustic, which is one of the two sets of caustics of 
the star-planet lens system located close to the primary star. 
Based on a similar perturbative approach, Asada (2002) pro- 
vides analytic expre ssions for the locations of the lensing im- 
ages. Bozza (2000) derived analytic expressions for the lo- 



cations of not only the central caustic but also the planetary 
caustic, the other set of caustics, which are located away from 
the central star. However, there has been no analytic work 
on the detailed properties of the caustics such as t he location, 
size, a nd shape, except the very recent work of Chu ng et alJ 
J2005I) (hereafter paper I) on the central caustics. 

Following paper I, we conduct a comprehensive and ana- 
lytic analysis on the properties of the planetary caustics. Un- 
der the perturbative approximation, we derive analytic expres- 
sions for the location, size, and shape of the planetary caustics 
as a function of s and q. Based on these expressions com- 
bined with those for the central caustics derived in paper I, 
we compare the similarities and differences between the plan- 
etary and central caustics. We provide an expression for the 
size ratio between the two types of caustics. We also derive an 
expression for the condition of the merging of the two types 
of caustics. We finally discuss the validity of the perturbative 
approximation. 

2. EMPIRICAL KNOWN PROPERTIES 

A planetary lensing is described by the formalism of a bi- 
nary lens with a very low-mass companion. Because of the 
very small mass ratio, planetary lensing behavior is well de- 
scribed by that of a single lens of the primary star for most 
of the event duration. However, a short-duration perturbation 
can occur when the source star passes the region around the 
caustics, which are important features of binary lensing. 

The caustics of binary lensing form a single or multiple 
closed figures where each of which is composed of concave 
curves (fold caustics) that meet at cusps. For a planetary case, 
there exist two sets of disconnected caustics. One 'central 
caustic' is located close to the host star. The other 'plane- 
tary caustic' is located away from the host star and its num- 
ber is one or two depending on whether the planet lies out- 
side (s > 1) or inside (s < 1) the Einstein ring. The size of 
the caustic, which is directly proportional to the planet detec- 
tion efficiency, is maximized when the planet is located in the 
'lensing zone', which represents the range of the star-planet 
separation of 0.6 < s < 1.6. The planetary caustic is always 
bigger than the central caustic. 
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3. ANALYTIC APPROACH 

We start from the formula of Bozza (2000) for the position 
of the planetary caustics (eqs. [49] and [50] of his paper). 
Keeping up to the first order term, the formula are expressed 
as 
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In these expressions, the coordinates are centered at the posi- 
tion on the star-planet axis with a separation vector from the 
position of the star of 
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where s is the position vector of the planet from the star nor- 
malized by 6*e (see Figs. [^and[2j. The origin of the coordi- 
nates corresponds to the center of the planetary caustic. For 
the pair of the planets with separations s and 1 /s, the cen- 
ters of the caustics are separated from the star by the same 
distance (because |r(s)| = |r(l /s)\) but directed toward oppo- 
site directions (because sign[r(s)] / sign[r(l/s)]). Therefore, 
the center of the caustic is located on the same and opposite 
sides of the planet with respect to the position of the star for 
the planets with s > 1 and s < 1, respectively. If one defines 
the lensing zone as the range of the planetary separation for 
which the planetary caustic is located within the Einstein ring, 
the exact range of the lensing zone is 
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To the first-order approximation, the size of the planetary 
caustic is proportional to q l l 2 as shown in equations Q and 
(|2ji. We will discuss the deviation of the approximation from 
the exact value in § 4.3. 

3.1. For Planets with s > 1 

In this case, between the two values of k in equation Q 
only the one with '+' sign is valid because the other one with 
'-' sign results in k 2 < 0. As a result, there exists only a single 
set of caustics for planets with s > 1 as shown in Figure^ 

The planetary caustic of the planet with s > 1 is composed 
of four cusps, with two of them are located on the £ axis and 
the other two are located on the r\ axis (see Fig.[0. The posi- 
tions of the individual cusps, (£ c , r) c )i, corresponds to the cases 
of sin 9 = (for the two cusps on the £ axis) and cos 9 = (for 
the other two cusps on the rj axis). Then, the positions of the 
cusps on the £ and r\ axes are expressed respectively as 
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FIG. 1 . — The location and shape of the planetary caustic of a planetary 
system with a star-planet separation greater than the Einstein ring radius (s > 
1 ). The upper panel shows the location of the star, the planet, and the resulting 
location of the planetary caustic. The coordinates are centered at the center of 
the planetary caustic, which is located on the star-planet axis with a separation 
vector r = s(l — 1/s 1 ) from the position of the star, where s is the position 
vector of the planet from the host star. The circle centered at the position of 
the star is the Einstein ring. The lower panel shows a blow-up of the region 
around the caustic enclosed by a box. Also marked are the definitions of the 
horizontal (Age) and vertical (Ar; c ) widths of the caustic and the designations 
of the individual cusps of the caustic. 



If we define the horizontal and vertical widths of the plane- 
tary caustic as the separations between the cusps on the indi- 
vidual axes (see Fig.0, the widths are expressed respectively 
as 
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where the expressions after the arrow are those evaluated to 
the first non-vanishing order in s in the limiting case of s ^S> 1 . 
Then, the vertical/horizontal width ratio is expressed as 
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In the limiting case of s ^> 1, A£ ~ A77 oc s~ 2 and 1Z C ~ 1, i.e. 
the caustic size decreases as s~ 2 and the shape becomes less 
elongated as the star-planet separation increases. 

3.2. For Planets with s < 1 

In this case, k in equation Q is valid only in the following 
range of 9 
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For 8 within these ranges, there are two possible values of k 
corresponding to the signs. As a result, there exist two sets of 



HAN & AN 



3 



0.4 
0.2 
PT 
0.2 
-0.4 



H~ -I 


\ 




- planet star 


1- 


H 



0.5 




FIG. 2. — The location and shape of the planetary caustic of a planetary 
system with a star-planet separation less than the Einstein ring radius (s < 1). 
Notations are same as in Fig. HI 



caustics for planets with s < 1; one above and the other below 
the star-planet axis (see Fig.|2j. 

Each of the caustics for the planet with s < 1 is composed of 
three cusps. One of them is located on the r\ axis but the other 
two are not located on either of the axes. The caustic meets the 
r) axis at T) c>+ = 2q l / 2 /[s(l +S 2 ) 1 / 2 ] and = 2q x l\\-s 2 ) x l 2 j s 
when cosf? = (see Fig.|2j. Among these two positions, the 
former corresponds to the cusp, and thus the location of the 
on-axis cusp is 

2q 1 ' 2 



0,±- 



sVT+ s z 

where the sign '±' is for the cusps located above and below 
the star-planet axis, respectively. 

If we define the vertical width of the caustic as the separa- 
tion between the two crossing points at r) Ct+ and rj c -, the width 
is expressed as 
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where the factor '1/2' is included into consideration that there 
exist two planetary caustics for planets with s < 1 and the 
expression after the arrow is that evaluated to the first non- 
vanishing order in s in the case of s -C 1. By defining the 
center of each caustic as the midpoint between the two cross- 
ing points (see Fig.|2j, its position is expressed as 
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The other two cusps occurs when d^/d9 = (or dr\ /d6 = 
0). This condition is satisfied when cos 2 2f? = l-3s 4 /4 (or 
sin 29 = 3s 4 / 4). Then, combined with the possible range in 
equation il Q , the values of 6 corresponding to the off-axis 
cusps are found to be 
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FIG. 3. — Upper panel: example planetary caustics of several planetary 
systems with different values of the star-planet separation .v and the planet/star 
mass ratio q. The orange and purple curves represent the loci of the center 
of the caustic as a function of s for planets with q = 5 X 10~ 3 and q = 5 X 
10 -4 , respectively. The coordinates are centered at the position of the planet- 
hosting star. The circle is the Einstein ring. Lower panel: separation between 
the caustic center and the position of the planet-hosting star as a function of 



With this value combined with equations Q and (|2j, the po- 
sitions of the off-axis cusps are expressed as 
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where kq = k(0q). In the limiting case of s <C 1, equation d!6t 
is approximated as 
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because ko — > s(\ +s 2 /4), sin^o 1, and cos6>o — > v / 3^ 2 /4. 
By defining the horizontal width as the separation between the 
two off-axis cusps, the width is expressed as 
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Once again, the factor '1/2' is included into consideration that 
there are two planetary caustics. From equations Jl 3I > and 
dl (St . the vertical/horizontal width ratio is expressed as 



d+s 2 r 1/2 -(i-^ 2 
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In the limiting case of s <C 1, each caustic shrinks as oc s 3 , c.f. 
oc s~ 2 for planets with s > 1, and K c ~ 4/3\/3 ~ 0.770, c.f. 
TZ C ~ 1 for planets with .s- > 1 . 

4. CAUSTIC PROPERTIES IN PLANETARY LENSING 
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Based on the analytic expressions derived in the previous 
section, we now investigate how the properties of the plan- 
etary caustics such as the location, size, and shape vary de- 
pending on s and q. We also compare the properties of the 
planetary caustics with those of the central caustics. 

4.1. Properties of Planetary Caustics 

In the upper panel of Figure [3] we present example plane- 
tary caustics of several planetary systems with different s and 
q. In the lower panel, we present the separation of the caustic 
from the planet-hosting star as a function of s. In Figure|4] we 
also present the variation of the caustic size (as measured by 
the horizontal and vertical widths) and the shape (as measured 
by the vertical/horizontal width ratio) as a function of s. 

The properties of the planetary caustics found from the fig- 
ures and the dependence of these properties on the planet pa- 
rameters are as follows. 

1. For s > 1, the location of the caustic center depends 
on s but not on q. On the other hand, for planets with 
s < 1, the caustic location depends on both s and q. In 
this case, the caustic is located farther away from the 
star-planet axis as q increases (see eq. 1141 ). 

2. Although the caustic size depends on the mass ratio as 
oc q 1 / 2 , the shape of the caustic does not depend on q 
and solely dependent on s (see eqs. 1101 and 1191 ). 

3. The rate of decrease of the caustic size with the increase 
of | log s | are different for planets with s > 1 and s < 1. 
Compared to the caustic of the planet with s > 1, the 
rate of decrease is steeper for the planet with s < 1. In 
the limiting cases of s ^S> 1 and s<Cl, the caustic sizes 
decrease as oc s~ 2 and oc s 3 for planets with s > 1 and 
s < 1, respectively (see eqs. 0, |9), \\3\, and fT8l ). 

4.2. Comparison with Central Caustics 



Chung et al. (2005) presented the analytic expressions for 
the location, cusp positions, width, and shape of the central 
caustics analogous to those presented in the previous section 
for the planetary caustics. The expressions for the location of 
the central caustic, analogous to equations Q and for the 
planetary caustic, are 



s+ 1 /s + 2(cos 3 </>-2cos0) 
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(20) 



(21) 



where (f> is a variable and the coordinates are centered at the 
position of the host star. There exists a single central caustic 
regardless of s and it has an elongated asteroid shape with four 
cusps, of which two are located on the £ axis and the other two 
are off the axis. The analytic expressions for the positions of 
the individual cusps, which are analogous to equations (|6} and 
for the planetary caustic with s > 1 and to equations d!2t 
and d!6i for the planetary caustic with s < 1, are 
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FIG. 4. — Variation of the size (normalized by q 1 / 2 ) and shape (as measured 
by the vertical/horizontal width ratio) of the planetary caustic as a function of 
the star-planet separation s. 

where cos0 t . = (3/4)(j+ l/s){l - [1 -(32/9)(s + 1/Vr 2 ] 1/2 }. 
The horizontal and vertical widths of the central caustic de- 
fined as the separations between the cusps on and off the star- 
planet axis are expressed respectively as 
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which are analogous to those in equations (|S} and (|5} for the 
planetary caustic with s > 1 and to equations dl 31 and dl St for 
the planetary caustic with s < 1 . Then, the width ratio of the 
central caustic is 

(s-l/s) 2 |sinV| 



(26) 



(s + 1 /s - 2 cos (j> c ) 2 ' 

which is analogous to those in equations dlOt and d!9l for 
the planetary caustics with s > 1 and s < 1, respectively. In 
the limiting cases of s ^ 1 and s<l, the size of the central 
caustic decreases respectively as 



{Aq/s 2 fors>l, 
\4^ 2 fors<l. 



(27) 



The planetary and central caustics have the following simi- 
larities and differences. 

1. Unlike the planetary caustic, the pair of the central 
caustics with separations s and 1 /s are identical as 
demonstrated by the fact that the inversion s 1/s in 
equations ( I20t and d2 II results in the same expressions. 

2. While the dependence of the size of the planetary caus- 
tic on the planet/star mass ratio is oc q 1 / 2 , the depen- 
dence of the central caustic is oc q. Therefore, the plan- 
etary caustic shrinks much more slowly with the de- 
crease of the planet mass than the central caustic. 
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between the two caustics is smaller than the half of the sum 
of the individual caustic widths, i.e. 



FIG. 5. — The size ratio between the planetary and central caustics as a 
function of the star-planet separation s and planet/star mass ratio q. For a 
representative quantity of the caustic size, we use the horizontal width. The 
region enclosed by the thick dashed lines represents the area in which the 
planetary and central caustics merge together, resulting in gradual obliteration 
of the distinction between the two types of caustics. 



3. For planets with s > 1, the rate of decrease of the size of 
the central caustic with the increase of | log s\ is similar 
to that of the planetary caustic with s > 1, i.e. A£ oc 
s~ 2 (see eqs. |8) and 1271 ). but smaller than that of the 
planetary caustic with s < 1, which shrinks as oc s 3 (see 
eq. ESI)- 

Then, what is the size ratio between the planetary and cen- 
tral caustics. If we use the horizontal width as a representative 
quantity for the caustic size, the size ratio between the two 
types of the caustics is found from equations (|SJl, dl 81 . and 
J24t and expressed as 
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fori > 1, 
-2 )cos6>o] for.? < 1, 
(28) 

where the additional subscripts 'p' and 'c' denote the plane- 
tary and central caustics, respectively. In Figure|5] we present 
the size ratio as a function of s and q. Since A£ c . c oc q while 
A£ C;P oc q 1 ' 2 , the dependence of the size ratio on the mass ra- 
tio is A£ CiC /A£ Cj p oc q x l 2 . For a given mass ratio, the size ratio 
is maximized at around s ~ 1 and decreases rapidly with the 
increase of | logs). 1 

As s — > 1, the location of the planetary caustic, i.e. r = 
s( 1 — 1 /s 2 ), approaches the position of the central star, around 
which the central caustic is located. Then the two types of the 
caustics eventually merge together, resulting in gradual loss 
of distinction between the two types of caustics. The condi- 
tion for the merging of the two caustics is that the separation 
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By using the analytic expressions for A£ c p (eqs. (3) and 1181 ) 
and A£ c c (eq. 1241 ). we compute the region of the caustic 
merging in the parameter space of (s, q) and presented in Fig- 
ure|5](the region enclosed by thick dashed lines). The region 
is confined in a small region around |s| ~ 1, but the width of 
the region increases as q increases because the caustic size 
increases with the increase of q. 

4.3. Validity of the Approximation 

Are the presented analytic expressions based on perturba- 
tive approximation good enough for the description of the 
caustics in planetary microlensing? We answer this question 
by comparing the two sets of caustics constructed based on 
analytic and numerical computations. 

In Figure |6] we present some pairs of the planetary caus- 
tics with different values of the planet parameters s and q. In 
each panel of the figure, the blue caustic is drawn by using 
the analytic expressions while the red caustic is the exact one 
based on numerical computations. For reference, we note that 
the mass ratios of the planets with masses equivalent to the 
Jupiter, Saturn, Neptune, and Earth around a host star with 
~ 0.3 Mq of the most probable Galactic lensing event are 
^~3x 10" 3 , 10" 3 , 2 x 10" 5 , and 10" 5 , respectively. From the 
figure, we find that although the deviation increases with the 
increase of the planet/star mass ratio, the analytic approxima- 
tion well describes the planetary caustic in most mass regime 
of planets (q < [0]1O~ 3 ). For the Earth-mass planet, we find 
that the two caustics are eventually indistinguishable. 

5. CONCLUSION 

We derived analytic expressions for the location, size, and 
shape of the planetary caustic as a function of the star-planet 
separation and the planet/star mass ratio under perturbative 
approximation. Based on these expressions, we conducted 
comprehensive analysis on the properties of the planetary 
caustics. Combined with the analogous expressions for the 
central caustics derived in paper I, we compared the similar- 
ities and differences between the planetary and central caus- 
tics. We also presented the expressions for the size ratio be- 
tween the two types of caustics and for the condition of the 
merging of the two types of caustics. These analytic expres- 
sions will be useful in understanding the dependence of the 
planetary lensing behavior on the planet parameters and thus 
in interpreting the planetary lensing signals. 
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FIG. 6. — Comparison of the planetary caustics based on the analytic (blue caustic) numerical (red caustic) computations for various values of the star-planet 
separation s and the planet/star mass ratio q. The coordinates are centered at the center of the individual caustics. The scales of the individual panels are set so 
that the caustics with the same s appear to have the same size. 



1 For the case of s < 1 , the change rate of the size ratio is reversed as | log s 
further increases beyond a critical value (| log .v ~ -0.3 or * ~ 0.5). However, 



this reversal occurs at the separation beyond the lensing zone. 
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